Spontaneous Symmetry Breaking in a 
Non-Conserving Two-Species Driven Model 



E. Levine f and R. D. Willmann | 

f Department of Physics of Complex Systems, Wcizmann Institute of Science, 
Rehovot, Israel 76100. 

| Institut fur Festkorperforschung, Forschungszentrum Jiilich, 52425 Jiilich, Germany 

Abstract. A two species particle model on an open chain with dynamics which is 
non-conserving in the bulk is introduced. The dynamical rules which define the model 
obey a symmetry between the two species. The model exhibits a rich behavior which 
includes spontaneous symmetry breaking and localized shocks. The phase diagram 
in several regions of parameter space is calculated within mean-field approximation, 
and compared with Monte-Carlo simulations. In the limit where fluctuations in the 
number of particles in the system are taken to zero, an exact solution is obtained. We 
present and analyze a physical picture which serves to explain the different phases of 
the model. 
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1. Introduction 

One dimensional driven diffusive systems have been shown in the last decade to 
demonstrate a variety of non-trivial types of behavior. In contrast to equilibrium 
systems with local interactions, driven systems were shown to exhibit boundary induced 
phase transitions, phase separation and spontaneous symmetry breaking even when the 
dynamics is local PUEJEJIIJE]- The latter phenomenon was first demonstrated for a 
driven two-species model with open boundary conditions, which became known as the 
'Bridge model' (HE]- This model is defined on an open lattice, where two species 
of particles, which can be thought of as carrying opposite 'charges', are moving in 
opposite directions. Although the dynamical rules which define the model are symmetric 
under charge exchange and left-right reflection, a mean-field approach and Monte Carlo 
simulations showed the existence of two phases in which this symmetry is broken. This 
result was subsequently supported by an exact solution for the limit of small extraction 
rates at the boundaries [Zj. 

Until recently, studies of one dimensional driven systems focused on conserving 
bulk dynamics. Non-conserving dynamics was only considered at certain sites, namely 
defect sites or boundaries. Recently, the existence of phase transitions in driven systems 



SSB in a non- conserving model 



2 



without bulk particle conservation was demonstrated, both for finite non-conserving 
rates ^U] and for rates which scale to zero with the system size [TIJEIJE]- 

In this work we study a variant of the Bridge model, in which non-conserving 
dynamics is introduced in the bulk of the system. All dynamical rates, both conserving 
and non-conserving, respect the CP-symmetry between the two species of particles. By 
its nature, this dynamics acts to balance the densities of the two species. If the non- 
conserving transitions occur at finite rates, spontaneous symmetry breaking does not 
occur. However, we find that when these rates are inversely proportional to the system 
size, spontaneous symmetry breaking appears. This model is related to the single- 
species non- conserving asymmetric exclusion process which was introduced and studied 
recently fTJEIJ^]. As in that model, the two-species model exhibits new phases with 
localized shocks. Transitions in the phase diagram can be understood by considering 
the position of the localized shocks. Furthermore, the phenomenon of induced localized 
shocks is observed, as predicted in [Tlj . 

The paper is organized as follows. The next section contains the definition of the 
model. In section |3] we calculate the phase diagram within mean-field approximation, 
for some cases where the bulk dynamics of the two species are decoupled. Different 
limits of the model are discussed, and in particular the phase diagram of the Bridge 
model is recovered as the non- conserving rates are reduced to zero. Next, we present in 
section 13] results of Monte-Carlo simulations, and compare them with the predictions of 
mean-field analysis. The phenomenon of induced shocks is addressed in section In 
section IH1 we present an exact solution for the case where fluctuations in the number of 
particles in the system are taken to zero. A physical picture is introduced and analyzed 
in section We conclude and summarize in section |H1 



2. Model Definition 



The model considered in this paper is defined on a one- dimensional lattice of size N. 
Each lattice site can either be occupied by a positive (+) particle, occupied by a negative 
(— ) particle, or vacant (0). The system evolves through three types of stochastic rates: 
In the bulk of the system particles move on the lattice according to the rates 

+ (Ao + +-^-+ 0-^-0. (1) 

In addition, each site in the bulk of the system, 1 < % < N, can change its state with 
rates 

+ +^0 -^0, (2) 

UI X UA U A 

corresponding to charge exchange, desorption of a particle from the lattice, and 
adsorption of a particle at an empty site. At the boundaries particles may be introduced 
and removed. At the left boundary, site i — 1, positive particles are introduced and 
negative particles are removed with rates 

0A+ -Ao (i=l), (3) 



SSB in a non- conserving model 



3 



while at the right boundary, i = N, negative particles are introduced and positive 
particles are removed with rates 

OA- +4o (i = JV). (4) 

Note that all dynamical rules, conserving and non- conserving, are CP-symmetric, 
namely symmetric under the exchange of positive-negative charges and left-right 
directions. Generalizations of this model to the case were both types of particles can 
move in both directions, and when the dynamical rules break the CP symmetry, will be 
considered elsewhere. 

Considering the bulk non-conserving rates, one distinguishes three possible 
scenarios [T^j. If the rates are finite, in the thermodynamic limit the bulk densities 
are dominated by the non conserving kinetics. Otherwise, if the rates decay to zero 
faster than 1/N, the bulk non-conservation should have no effect, and the properties of 
the system in the thermodynamic limit are identical to those with bulk conservation. 
The third case, which we consider here, is the one in which the non- conserving rates scale 
down linearly with the system size. It is useful to introduce the notation uja = £l/N, 

Without the non-conserving dynamics in the bulk of the system, Eq. (j2J, the model 
is identical to the Bridge model introduced in [B] and further studied in jHIHlE]- Let us 
summarize first the known results for this model. At q = 1 mean-field analysis predicts 
the following phase diagram !§.]. Two phases which obey the CP-symmetry of the model 
are identified: a maximal-current phase, where the bulk density of each species equals 
|, and a low-density phase, where the bulk densities of both species are identical, both 
smaller than |. In addition, the model exhibits two phases in which the symmetries are 
broken. In one phase, termed the Low-High phase, one of the species, spontaneously 
chosen, sustains a bulk density larger than ~, while the bulk density of the other is 
smaller than |. In the other phase both species sustain low albeit non-equal densities. 
This phase spans a very small area of the mean-field phase diagram, and is very difficult 
to recover in Monte-Carlo simulations. Consequently, the existence of this low-density 
asymmetric phase under noisy dynamics has been the subject of some debate |H1IE]- 

The model considered here can be thought of as two interacting single-species totally 
asymmetric exclusion processes with bulk non-conservation (NC-TASEP). The NC- 
TASEP is an extension of the well-known totally asymmetric exclusion process (TASEP) 
[TH] with bulk absorption and desorption dynamics. Using mean-field calculations [12], 
which were argued to be exact [T3|, the phase diagram of this model was obtained. In 
addition to the three phases of the TASEP, namely the maximal-current phase, the 
high-density phase and the low-density phase, it was found that the NC-TASEP may 
also exhibit four additional phases. The most interesting one is a Shock phase which 
consists of a localized shock in the bulk of the system, separating a low-density region 
to its left from a high-density region to its right. In the bulk-conserving TASEP shocks 
appear only on the boundary line between the high-density and low-density phases. On 
that line a delocalized shock appears in the system. As the position of the shock is 
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equally probable at any site in the system, the average profile on the transition line is 
linear. In contrast, the NC-TASEP exhibits a distinct phase in which a localized shock, 
whose position is selected by the dynamics, is present. The existence of this phase plays 
a main role in our analysis of the two-species model. 

The maximal-current phase of the TASEP appears in the NC-TASEP only when 
the two non-conserving rates, namely particle absorption and desorption, are equal. In 
this case there exist three more phases. These include a Low-Max phase, in which the 
density increases linearly from a boundary density < \ towards density |, where it 
remains constant for the rest of the system ; a Max-High phase, in which the density 
rises linearly from a constant profile of density \ to a boundary density > \ ; and a 
Low-Max-High phase, in which the density rises linearly from a left-boundary density 
< | towards |, where it remains constant up to a point where it climbs linearly again 
towards a right-boundary density > \. 

3. Mean Field Theory 

In this section we study the mean-field equations of our model in the thermodynamic 
limit N — > oo. We introduce the occupation variables and 9i, such that Tj = 1 (9j = 1) 
if site % is occupied by a positive (negative) particle, and otherwise. The densities of 
the positive and negative particles are then defined by 



Pi = (n) mi = (9i) , 

where angular brackets denote averaging over realizations. 

The time evolution of the particle densities is governed by 



(5) 




(6) 



where the currents are given by 



3? = 

k = (9 i {l 




q)0i+i)) 

q)n-i)) 



(7) 



and the source terms are 



= uj a (1 - pi - rrii) - uj D pi + lo x (m - p t ) 
S~ = u A (1 - rrii - p^ - LU D irii + lo x (Pi - m) ■ 
At the boundaries the source terms vanish, and the currents are given by 



(8) 




Ppn 

f3m\ 



Pi - mi) 



(9) 



3n+i 



p N - m N ) . 
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The mean-field theory for this model is defined by replacing two-point correlation 
functions with products of one-point averages. Within this approximation, the currents 
become 

3i =Pi( 1 ~ Pi+i ~ (! - q) m i+i) 

j7 = m (i - m-i - (i - q)Pi-i) ■ (io) 

In the steady-state all time derivatives vanish, and one has 

3t=3ti + Sf jr=3t+i + St- (11) 

Defining J t + = jf - Y!k=o S k $i = 37 ~ J2k=i S k> one notices that in fact J t + = J + 
and = J~ are conserved throughout the lattice, and J + = Jq, J~ = j'^r+i- 

3.1. Solution of the mean-field equations in the bulk- decoupled case 

The case q = l,ux = 1 is special, as in this case the bulk equations (JHJ) with the 
mean-field currents are decoupled. The currents and source terms in this case are 
just 



37 = Pi^-Pi+i) S+ = %[l-(2 + u D ) Pi ] 

jr = mi{l- mi_i) = § [1 - (2 + u D ) mi 



(12) 



In the bulk of the system, the hopping rates for say a positive particle do not distinguish 
between a negative particle and a vacancy. Also, the fact that attachment of a positive 
particle is limited by the presence of negative ones, is exactly compensated by the charge 
exchange rate. The coupling between the two species is limited in this case only to the 
boundaries. Following [B], one readily notices that upon the definition 

a(l -px- mi) _ J + 3o 



a: 



i-px ^ + ^F # + ^ 



(3 



_ = a(l-p N -m N ) = J Jq n „v 

a l-m N J^ + J±±s± i^.ii' 1 } 

a P a + 

with S ± = Y^hLx the problem is reduced to two single-species NC-TASEP. One 
process corresponds to the positive particles with injection rate a + at the left boundary 
and ejection rate (3 at the right, and the other to the negative particles with injection 
rate a~ at the right boundary and ejection rate (3 at the left. The two processes may 
or may not share the same phase. The latter case corresponds to a phase of the two- 
species system, where the symmetry between the two species is broken. In the other 
case, it may be that the average densities of the two species are not equal, although 
the two lie in the same phase of the corresponding NC-TASEP. A trivial restriction on 
the possible phases in the model is pi + < 1 at all sites. This immediately excludes 
several possibilities, such as ones which mix the high-density phase of one species with 
anything but the low-density of the other. 

In this section we explore the possible phases in the bulk decoupled case. Symmetric 
phases are presented first, followed by asymmetric phases. In the symmetric phases, 
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Figure 1. Density profiles of the symmetric phases, as obtained by integrating the 
mean-field equations for a system of size N — 128. (a) Max Phase, a — 3.0, (3 = 0.8, 
= 0.2, q = 1. (b) Low-Max Phase, a = 1.0, (3 = 0.7, fi = 0.2, q = 1. (c) Low 
Phase, a = 0.1, (3 = 0.8, n = 0.2, q=l. (d) Low-Max-High Phase, a = 5.0, = 1/3, 
ft = 0.5, g=l. (e) Shock Phase, a = 3.0, /3 = 0.25, ft = 0.2, q = l. 



a + = a~ , so for these phases only a + is quoted in the following. For the asymmetric 
phases it is always assumed, with no loss of generality, that the positive particles are in 
the majority When describing density profiles we always take a language in which the 
lattice is rescaled to the segment [0, 1]. The emerging phase diagram is discussed in the 
following subsection. 

3.1.1. Maximal- current symmetric phase. In this phase the bulk density of both species 
is |, and the boundary currents are given by 

Jo = Jo = \ ■ (14) 
The conditions for this phase to exist are 

« + > ~ P>\. (15) 

Typical density profiles of the two species in this phase, as obtained from integrating 
the mean- field equations, are shown in figure ^ (a). 

3.1.2. Low - Max symmetric phase. The density profile in this phase is composed of 
a low density part where the density increases linearly with slope Q on the rescaled 
lattice, as well as a part with constant density ~. The boundary currents are 

j + = a+(l-a + ) j ~ = 7- (16) 
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The conditions for the existence of this phase are 

a + < 1 - a + > l --tl /?<-. (17) 

The density profiles shown in figure ^ (b) for a finite system furthermore exhibits a 
boundary layer, which does not scale with the system size. 

3.1.3. Low density symmetric phase. In this phase both densities remain below |, 
increasing throughout the system with constant slope Q (figure life)). The boundary 
currents are given by 

j+ = a+(l - a + ) jo = (a + + 0)(1 - a + - O) . (18) 

Necessary conditions for the existence of this phase are 

a+ < 13 - Q for (3 < ~ 

a + < ^ - ^ for V > g • ( 19 ) 

Inserting the boundary currents (fTHj) into (JT3Jl yields a quadratic equation for a + . Using 
fTT)j) one readily identifies the relevant solution. 

3.1.4- Low - Max - High symmetric phase. The density profiles in this phase are a 
mixture of three different pieces - a linear profile of low densities and constant slope Q, 
followed by a flat density profile at density |, and a linear profile of high-densities of 
the same slope (figure dd)). The boundary currents are now 

j+ = a + (l-a + ) 3 - =f 3{i-(3). (20) 

This phase region is defined by the conditions 

a + > l --n &<\. (21) 



Again, one solves the equation for cr given by |T3j) and (JzOJ) . and uses (|2T|) to identify 
the relevant solution. For a finite system (as seen in figure dd)) the transition between 
the three parts is not sharp. It only becomes so on the rescaled lattice as iV — >• oo. 

3.1.5. Shock symmetric phase. This phase is characterized by a localized shock, 
separating a low-density region from a high-density one, both of linear profile with 
slope Q. The boundary currents are given by 

j+ = a + (l-a + ) j-=(3{l-(3). (22) 

Notice that a + in this phase is identical to the one of the Low-Max-High phase. The 
Shock symmetric phase is defined by the conditions 

(3-Q<a + <^-Q. (23) 
The position of the shock x s is given by 

- , I _ 2/3 -(! + «) + y/(l + a) 2 - 4a/3 1 

Xs ~ 2il + 2 ~ 4fi + 2 ' 1 j 
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Figure 2. Density profiles of the asymmetric phases, as obtained by integrating the 
mean-field equations for a system of size N = 128. (a) Shock-Low asymmetric Phase, 
a = 0.1, (3 = 0.05, n = 0.02, q=l. (b) High-Low asymmetric Phase, a = 1.5, (3 = 0.1, 
tt = 0.02, q = l. (c) Low asymmetric Phase, a = 0.82, (3 = 0.31, ft = 0.02, q = 1. 



In contrast to the thermodynamical limit, the shocks in a finite system such as in figure 
[U (e) are not sharp. 

3.1.6. Shock - Low asymmetric phase In this phase the majority species exhibits a 
localized shock, while the minority species is in the Low phase throughout the system 
(figure El (a)). The boundary currents for the two species are 

j+ = a+(l - a + ) Jo = + - aT - SI) 

j+ = P(l-P) S = a-(l-a-). (25) 

The conditions for the existence of this phase are given by 

a + > p-n a + <p + n (3< \- ( 26 ) 

The equation for a~ does not involve a + , and can be solved as in previous phases. 
Plugging this solution into the equation for a + one can solve the equation, and identify 
the only solution which obeys (j2fij) . 

3.1.7. High - Low asymmetric phase This phase is analogous to the strong asymmetric 
phase of the Bridge model. In this phase the majority sustains a high density in the bulk 
of the system, while the minority density is low. Here, however, the density profiles are 
not constant, but rather of opposite slopes ±Q (figure H| (b)). The boundary currents 
for the two species are 

j + = 03 + 0) (1 - P - fi) jo = (a~ + 0)(1 - a' - ft) 

j+ = P(l-P) j N = a -(l-a-). (27) 

The conditions for the existence of this phase are 

a + >P + Q p + Q<^. (28) 

Expressions for a 1 * 1 are obtained from ()27|). PH|) just as in the Shock-Low phase. 
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3.1.8. Low asymmetric phase In this phase both particle species maintain a low density 
profile with constant slope Q. Still, the phase is asymmetric as the boundary densities 
of the two phases are different. An analogous phase is also observed on the mean field 
level in the Bridge model. The boundary currents are given by 

j + = a + (l - a+) j - = (a- + - a - fi) , 

j+ = a"(l - a") j+ = {a + + - a+ - O) . (29) 

Plugging the currents (j2El) into (fHfy gives 

a + (l-a + ) (a~ + - a~ - n) 



a + = 1 - 



a (3 



a-(l-a-) ( ct + + ^)(l- a + -^) 
a =1 (30) 

a p 

Let D = a + — a~ and 5 = a + + a~ . Subtracting the equations in ([3~U|) yields 

D = D ( {1 - S) ^™). (31 ) 



a/? /? 

In the asymmetric phase, D ^ 0, thus an expression for S is obtained. Summing the 
equations in (jHUj) and using this result yields D as a function of a and (3. The effective 
boundary rates are obtained as a + = ^(S+D) and a~ = ^(S—D). Necessary conditions 
for the existence of this phase are 

a + < (3 - Q and (3 < - 

a + < - — Q, and > - 
2 H ~ 2 

D > . (32) 
A typical profile for this phase is shown in figure 121 (c). 

3.2. Phase Diagram 

In the previous section we have listed all phases which are found to exist in this model, 
and derived the phase boundaries in which they reside in parameter space. We now 
turn to describe the emerging phase diagram. First note that the full parameter space 
is covered by the four symmetric phases. In fact, all asymmetric phases reside in regions 
of phase space where the low-density symmetric state is also stationary. Which of the 
solutions is realized is a matter of stability, as will be discussed shortly. Except for the 
boundary of the Maximal- current phase, the phase boundaries of the symmetric phases 
share a common point Q in the (a, (3) plane, given by 

Q-( l -^,l)- (33) 



4Q 2 y 

The intersection point of the phase boundary between the Low and the Shock symmetric 
phases meets the (3 axis at the point 

R=(0,Q) . (34) 
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Figure 3. Mean-field phase diagram for the bulk-decoupled case with lod = 
and f2 = 0.2 (left), = 0.02 (right). The inset focuses on the regime where the 
low-asymmetric phase is most pronounced. Abbreviations used for the phases: M - 
maximal-current symmetric phase; LM - Low-Max symmetric phase; LMH - Low-Max- 
High symmetric phase; S - Shock symmetric phase; L - Low symmetric phase; SL - 
Shock-Low asymmetric phase; HL - High-Low asymmetric phase; LL - Low asymmetric 
phase. 

In figure|2]we plot the phase diagram for the cases Q = 0.2, 0.02. Note that at the higher 
value of Q no asymmetric phases exist, while at the smaller value all phases described 
in the previous section exist. This fact will be addressed below. Taking Q to zero the 
original phase diagram of the Bridge model is retrieved, as discussed in subsection 13. 2. HI 

3.2.1. Stability. As mentioned above, all asymmetric phases reside in regions of phase 
space in which the symmetric low density is also a stationary solution of the mean-field 
equations. On the mean field level, the realization of one stationary solution rather 
than the other is a matter of initial conditions. In all cases both the symmetric and 
asymmetric solutions are linearly stable. However, any initial condition for which the 
density of at least one of the two species is higher than | in some region evolves into 
the state of broken symmetry. Thus a disordered initial condition, in which the density 
of particles at any site is an independent uniformly distributed random variable, resides 
in the basin of attraction of the asymmetric solution. 

Considering the model beyond mean-field approximation, where the dynamics is 
noisy, one expects a random perturbation to take the system away from the symmetric 
solution. In physical terms this can be understood by the fact that nucleation of a high 
density domain leads to its flow towards the boundary, where it reduces the inflow of 
particles of the other species due to the exclusion interaction. Once the symmetry is 
broken, the high density phase expands and takes its steady state position. This picture 
is substantiated in a quantitative manner for the limit B — > in section El Note that 
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this line of argument cannot be followed for the asymmetric Low phase. 



3.2.2. Extinction of the asymmetric phases at high Vt. So far we have considered the 
case ujd = 0, where detachment of particles from the bulk is suppressed. In this case 
the non-conserving rates allow for attachment of particles of either species with rate uoa 
and for charge exchange with rate uox- While the former process affects the densities of 
both species in the same way, the charge exchange process tends to diminish the density 
difference between species. Thus it is clear that for large ujx ~ fi, when this process 
becomes dominant, the asymmetric phases will vanish. This effect can only increase 
in the presence of detachment, which acts stronger on the majority species. In the 
case up = we find that the High-Low phase ceases to exist beyond Q ~ 0.035. The 
Shock-Low phase vanishes at Q ~ 0.138, and with it the Low asymmetric phase. The 
vanishing of the asymmetric phases can be understood in a more quantitative manner 
from the blockage picture described in section 



3.2.3. The limit — > 0. In the limit of f2 — > the non-conserving model considered 

here must coincide with the Bridge model. In this limit, the point R defined above is 

shifted towards the origin. The point Q is pushed to infinity, which means that the 

Low-Max-High symmetric phase cannot exist. The boundary line between the Low 

symmetric and the shock symmetric phase is given generally by 

i i , 1 _ AO 2 

Pls = n + - - -Vl-4an for a < — ^— . (35) 

Thus /?ls — * for all a as Q — ► 0, and the shock symmetric phase cannot exist. 
Furthermore the Low-Max symmetric phase vanishes in this limit. The boundary 
between this phase and the Maximal- current phase, 

coincides in the limit SI — ^ with its boundary with the Low phase, 

a r 1 l-4fi 2 
P = ^7 T-T- for - < a < — . (37) 

Thus, of the symmetric phases only the Maximal- current phase and the Low phase 
remain in the Q —>■ limit, as expected from the Bridge model. 

As for the asymmetric phases, one notices in the same way that as Q decreases, 
the High-Low phase region grows on the expense of the Shock-Low phase. As Q — > 0, 
the boundary lines of both phases coincide, and the Shock-Low phase ceases to exist. 
Exact expressions for these phase boundaries are rather lengthy, and are omitted here. 
Finally, the low density asymmetric phase takes in the Q — > limit its form as in the 
Bridge model. 
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Figure 4. Mean-field phase diagram for the bulk-decoupled case, with £1 = 0.02 and 
Up = 1 (gray lines) , ujj = 2 (dark lines) . The boundary line between the asymmetric 
low-density phase and the symmetric Low phase is omitted. Abbreviations for the 
different phases are the same as in figure 



3.3. Detachment from the bulk: the case iod ^ 

We now turn to consider a more general case, still within the bulk-decoupled regime 
(q = l,ou D = ujx), where the non-conserving dynamics includes detachment of particles. 
This case corresponds to the NC-TASEP with non-equal attachment and detachment 
rates. The phase diagram of that model includes only three phases: high density, low 
density and localized shock. All phases where a part of the profile is constant at ~ do 
not exist. This comes from the fact that the equilibrium density of the kinetics is given 
by the Langmuir density rather than |. 

The phase diagram of our model in the case u>n ^ exhibits only two symmetric 
phases - the symmetric low-density phase and the shock symmetric one. All asymmetric 
phases which were obtained for the case ujd = are also present here. As in the NC- 
TASEP, the density profiles in this case are not linear, but rather curved. This makes 
the analysis of the phase diagram somewhat more cumbersome, although not different 
in principle from the one presented in section 13.11 Details of this calculation are given 
in the appendix. 

In figure |U we plot the mean- field phase diagram for the cases up = 1 and up = 2, 
where up is defined as before by «b = lob/^a- The boundary line between the 
asymmetric low-density phase and the symmetric Low phase is not presented. We 
could not obtain the boundary densities in this phase in a closed form. For several 
values of up we have found numerically that this line lies just above the transition line 
between the Shock-Low phase and the asymmetric Low phase. Certainly, the region of 
phase space covered by this phase does not increase compared with the case ojd = 0. 
It is readily noticed that the part of phase space spanned by the High-Low asymmetric 
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phase decreases as ujd is increased. This is expected from the fact that the detachment 
process acts stronger on the majority phase, thus reducing its density. For any given Q 
the detachment process can be increased to a level in which the High-Low asymmetric 
phase does not occur. 

The detachment process can be considered as cooperating with the boundary 
ejection rate (3, and competing with the boundary injection rate a. It is no surprise 
then that the asymmetric Shock-Low phase grows in the a direction of phase space 
and shrinks in the (3 direction as ujd is increased. In the regions of phase space which 
compose the three phases identified only for the case ojd = 0, the density profiles 
approach continuously, as is decreased towards zero, to the ones described in the 
corresponding phases of the ojd = case. 

4. Monte-Carlo simulations 

As was already noticed in earlier works about the Bridge model [313, the mean field 
phase diagram captures the correct topology of the phase boundary lines. The exact 
location of the boundary lines, however, is shifted in the noisy model with regard to 
those of the mean field solution. For our model we did not try to obtain the exact 
location of the phase boundary lines from Monte-Carlo (MC) simulations. Still we note, 
based on our simulations, that these lines cannot lie too far from those of the mean-field 
phase diagram obtained in the previous section. Here we concentrate on giving evidence 
for each of the phases by finding representative points in parameter space. In figure 03 
we present the density profiles of the two species in the different phases as obtained from 
MC simulations. The profiles in each phase were taken at the same parameters as the 
respective mean field profiles shown in figures Q and EJ For all phases shown here the 
mean field profiles capture the features of the noisy model. While MC simulations were 
done for a system of size N = 2000, the mean field results are obtained for N = 128. A 
quantitative comparison between the profiles becomes meaningful in the limit iV — > oo 
using rescaled coordinates x = i/N. In this limit the localized shocks in the Shock-Low 
and symmetric Shock phase become sharp jTT] . In figure El the density profile for the 
Low asymmetric Phase is omitted. This phase is addressed in the following subsection. 

4-1. The Low asymmetric phase 

The existence of the Low asymmetric phase is an issue of longstanding discussion jHllHl- 
It was noted already in the mean-field solution of the Bridge model jB] that the region in 
phase space covered by this phase is very small compared to the others. Furthermore MC 
simulations indicate that the particle densities in this region of phase space fluctuate 
strongly These facts also hold true in the model considered here. Therefore we 
refrain from presenting a MC density profile as was done for the other phases. In jH] the 
existence of the Low-asymmetric phase in the Bridge model is deduced from sampling 
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Figure 5. Density profiles as obtained from a Monte-Carlo simulation of a system of 
size N = 2000. (a) Max Phase, a = 3.0, (3 = 0.8, fl = 0.2, g = 1. (b) Low-Max Phase, 
a = 1.0, 13 = 0.7, n = 0.2, q=l. (c) Low Phase, a = 0.1, (3 = 0.8, ft = 0.2, q = 1. (d) 
Low-Max-High Phase, a = 5.0, = 1/3, tt = 0.5, g = 1. (e) Shock Phase, a = 3.0, 
13 = 0.25, ft = 0.2, q = 1. (f) Shock-Low asymmetric Phase, a = 0.1, /3 = 0.05, 
ft = 0.02, g = l. (g) High-Low asymmetric Phase, a = 1.5, (3 = 0.1, ft = 0.02, q=l. 



the probability distribution P(p, m), where p and m are the average densities of positive 
and negative particles in the system respectively. In our model, this line of argument 
fails. On the level of average densities, the Low-Low and Shock-Low phases cannot be 
distinguished. This is because both phases can exhibit distributions P(p, fh) with two 
peaks at p and m smaller than |. 

The blockage picture outlined in section [7| yields no arguments in favor of the Low 
asymmetric phase. At the upper boundary of the Shock-Low phase the localized shock 
position retracts to i s = 0. This allows particles of both species to enter the system. 
A symmetry broken Low phase beyond this point would require some kind of blockage 
being formed at the exit of the majority species. The nature of such a blockage is not 
clear. It remains to be clarified, whether the existence of this phase, which is evident in 
the mean field treatment, can be demonstrated in the noisy model. 
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Figure 6. Induced shocks. (Left) Density profiles in the symmetric shock phase, 
as obtained in Monte-Carlo simulations. A primary shock in one of the species is 
accompanied by an induced shock in the other species. Here q = 2 and a = 3.0, 
j3 = 0.25, £1 = 0.02. (Right) Density profiles in the Shock-Low phase, as obtained in 
Monte-Carlo simulations. The localized shock in the majority species is accompanied 
by an induced shock in the minority species. The profile of the minority species is 
additionally shown when multiplied by a factor of 30, in order to demonstrate the 
induced shock phenomenon. Here q = 2 and a — 0.1, j3 — 0.05, fl — 0.02. 



5. Induced shocks 

In the case q ^ 1 the bulk dynamics of the two particle species are not decoupled. Thus, 
one cannot solve the mean field equations in the way it was done in sectional Still, the 
phase structure can be explored by integrating the mean field equations numerically, or 
by MC simulations. It shall not be attempted here to give the full phase diagram of 
the model. We do note, however, that phases with broken symmetry exist also in the 
general case. 

In the Shock symmetric phase and the Shock-Low phase the coupling of the 
dynamics of the two species gives rise to an induced shock phenomenon. Here the 
existence of a localized shock in the density profile of one species induces a shock in 
the density profile of the other species. For example, in the Shock symmetric phase one 
notices that the density profile of each one of the species exhibits actually two shocks 
in the steady state (figure EJ). One is a primary shock, created by the same localization 
mechanism which was already identified in the bulk-decoupled case. The existence of 
this shock, albeit not its detailed properties, relies only on the properties of the density 
profile of the very same species. The second shock is induced by the primary shock in the 
density of the other species, and it shares its location. Across both shocks the current 
is continuous. This phenomenon also occurs in the Shock-Low phase (figure EJ). Here, 
the localized shock in the majority phase induces a shock in the Low phase. In fact, 
it was shown in ^3] that for general two species systems with coupled density-current 
relation a shock in one particle species induces a shock in the other one. 

For general q the current-density relation of, say, the positive particle species, 
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j^ipi,Tni), depends on the local density of both species. In general this current-density 
relation is not known, except for two cases: the case q — 1 (the decoupled case), where 
j + = p(l — p), and the case q = 2, where j + = p{l — p + m) 

As discussed in [T3] the current across a localized shock is continuous. This 
requirement implies when q = 1 that shocks are symmetric with respect to p = 1/2, 
irrespective of the local density m of particles of the opposite species. In general, 
however, this is not the case. The properties of both the primary and induced shock in 
the density profile of each species rely on the local densities of both. 

The continuity of the current across the shocks can in principle be used to determine 
the properties of the primary and induced shocks, if one can develop the density 
characteristics from the boundaries of the system. For the case q = 1 the equations 
for the two density profiles are decoupled, and one uses this method to determine the 
position of the shock. Of course, no induced shocks are present in this case. For general 
q, however, the equations for the density profiles are coupled, and an analytical solution 
is not available. 

6. Exact solution for the limit (3, — > 

In [7] a toy model was presented to allow for the exact solution of the limit (3 — ► of 
the Bridge model. In this section a generalization of the toy model is presented. The 
solution of this model gives an exact description of the (3 — > 0, Q — > limit of the model, 
and proves that to lowest order in f3 mean-field theory recovers the exact phase diagram. 

In the limit (3 — > 0, Q — > the only relevant configurations are those composed of 
three blocks, containing (from left to right) negative particles, vacancies, and positive 
particles. A configuration of this type is long-lived, as all exit rates from it scale to zero. 
In this limit, all other configurations rearrange themselves into one of these three-block 
configurations. A configuration of this type is identified by two variables, y and x, 
defined as the size of the left (negative particle) block and the right (positive particle) 
one, respectively. Thus, for example, 



Let us assume that the system is in a three-block configuration (x, y) , and consider 
the ways it can leave it. First, a particle can leave the system through a boundary 
with rate (3, leaving a vacancy behind it. This vacancy can travel into the system 
with rate (1 + a)" 1 , in which case the system is again in a three-block configuration. 
Otherwise, the particle which had left the system can be replaced by a particle of the 
opposite species with rate a(H-a) -1 . On a short time scale this particle travels through 
the system until it joins the block on the other side, thus returning the system into a 
three-block state. 

Another possible way out of a three-block configuration is through the non- 
conserving processes in the bulk. First, a particle can be attached to the system in 
the vacancy domain with rate ua{N — x — y). This particle joins on a short time scale 
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Figure 7. Toy model for the limit /3, 17 — > 0. (Left) The rates defining the 
corresponding random- walk model i|38|) . (Right) Flow fields of the model. 



to the block of its own species. Second, a positive (negative) particle can be detached 
from the system with rate oje>x (upy), thus creating a vacancy within a particle block. 
On a short time scale this vacancy travels into the system and joins the vacancy block. 
Finally, a positive (negative) particle can change its species with rate u x x (u x y), and 
move from one particle block to the other. 

When the last particle of its species leaves the system, the other type of particles 
can rush into the system through the boundary. The system fills rapidly with particles 
of this type. Thus, the only possible configurations with x = or y = are (0, N) and 
(iV,0). ' 

To summarize, consider a two-dimensional random walker, whose position (x, y) 
corresponds to the block-configuration of the two-species system. The transition rates 
for this walker are (see also figure EJ) 

with rate loa{N — x — y) 
with rate u>a(N — x — y) 



x,y) 
x,y) 
x,y) 



(x + I,?/) 
(x,y + 1) 

(x - l,y) 



with rate f3 



+ oj d x 



x,y) ->(x,y- 1) 
x,y) -> (x- l,y+ 1) 

x,y)-*{x + l,y- 1) 

l,y)^(0,N) 
0) 



with rate 
with rate 

with rate 

with rate 
with rate 



1 + a 

Pt-. 1- u D y 

1 + a 

(3— h ujdx 

1 + a 

n a 

P— h to D y 



(38) 



1 + a 



It can be shown that this toy model is mapped exactly to the two-species model in the 
limit (3, Q — > 0, , in the sense formulated in appendix A of [7]. 

Let us first consider the case where the dynamics defined in ()38|) leads to a fixed 
point solution. This is the case where the net flows on both the x and y directions 



SSB in a non- conserving model 



18 



vanish at some point (x*,y*). The fixed points must satisfy the equations 



(uj x - uj a )x* - (uj x + u A + u D )y* = - u A N + - 




(39) 



whose solution is 



^ _ if __ l A _ g \ 



N N 2 + u D V 0(1 + a) J 



(40) 



Interestingly, the fixed point does not depend on the charge exchange rate ujx- Notice 
that this fixed point can only exist if < x*/N,y*/N < I. Indeed, x* and y* of (JUjJ) 
always meet the second condition. The first condition, however, is only met for 



Otherwise, the random walker is always biased towards the axis of the (x, y) plane which 
is closer to its position (see figured). In this case, Q < (3 — > 0, one recovers the toy 
model of [Zj which yields a stable state of broken symmetry. One of the species then 
occupies most of the lattice, corresponding to the High-Low asymmetric phase of the 
model. 

Thus, the toy model yields two phases in the limit (3,Q —>■ 0. For (3 < Q(l + a) one 
has a symmetric phase, with a dominating three-block configuration described by (|40|) . 
This corresponds to a symmetric shock phase in the model, where the shock position of 
the positive particles is x s = 1 — x* /N, with x* given by (|4TJ|) . Otherwise the system is 
in the High-Low asymmetric phase, with the line (3 = 0(1 + a) serving as the transition 
line between the two phases. 

For the bulk-decoupled case, u A = ojx and q — 1, it is illuminating to compare these 
exact results with the ones obtained by mean field. The mean field analysis, performed 
in section El predict in the limit /3, Q — *■ the two phases obtained in the toy model. 
The mean field transition line between the two phases ()23)1 is identical, to first order in 
P, to the line (3 = Q(l + a) of the toy model. Also the shock position x s calculated in 
the toy model is identical to first order in (3 to the one (jUj) calculated in mean field. 
This result also holds in the case ujd ^ 0. We thus conclude that the mean- field solution 
is exact to first order in (3. 

7. Blockage picture 

In this section we combine the mean-field and stability analysis, the simulation results, 
and the toy model into a physical picture. Following 8\ we term it the blockage picture. 

Qualitatively, typical configurations in the asymmetric phases can be described 
in terms of blocks of the two species, which spread from the 'exit' boundary into the 
system (figure |Hl and compare the toy model description in the previous section). The 
density profile within each block is not constant, but this feature is not relevant here. 
As mentioned above, a block of one species stalls the entry of particles of the other 



p<n{i + a). 



(41) 
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Figure 8. Blockage picture. Schematic picture of the instantaneous density profiles. 
Density profiles of positive particles are depicted by dark lines, negative particles by 
gray lines. Here, as in the text, we assume that the positive particles are the majority 
species. 



species through the boundaries, thus serving as a blockage. The possibility of particles 
to enter the system in the bulk serves to stabilize the domain size. 

In the High-Low phase the block of the majority species covers the entire system, 
while in the Shock-Low phase the block is limited to some part of the system. The 
fluctuations in the size of this block, corresponding to the width of the localized shock, 
are limited to an area of size ~ iV -1 / 2 [11J. The minority block in both phases is 
unstable in the sense that the domain wall between it and the bulk region drifts towards 
the boundary. Averaging over the positions of the domain wall results in the exponential 
decay of the mean field density profile from the left boundary. 

In |Sj it was observed that an instantaneous configuration in the low-density 
asymmetric phases comprises a small block of the majority phase, which is limited 
to the vicinity of the boundary. The formation of this block prevents particles of the 
other species from entering the system, thus leading to symmetry breaking. However, 
this block does not survive for times which are exponential in the system size. 

The block picture is extended into the symmetric shock phase. Here the two blocks 
are covering equal distances from the 'exit' boundaries. The sizes of the two blocks are 
again macroscopic and localized, in the sense that the size fluctuations vanish as N~ l l 2 . 

We now turn to describe the different phases of the model in the language of block 
configurations. To this end we take a stroll along a line of constant a in the phase 
diagram, starting from the symmetric shock phase and going up in (3. This line is 
chosen such that it cuts through all asymmetric phases. 

In the shock symmetric phase, the two blocks inhibit, in a symmetric way, the inflow 
of particles of the opposite species. Increasing j3 decreases the sizes of the two blocks. 
However, as long as their size is macroscopic, the blocks keep their role of lessening the 
ability of particles of the other species to enter the system. 
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P P 

Figure 9. Blockage size as a function of 0, as calculated in mean-field for the bulk- 
decoupled case. Here a = 0.2, lud = and (a) = 0.02, (b) £1 = 0.03. 

As the boundary line is approached, the size of the blocks is reduced to zero. 
Now the road is open for both species to enter the system. Due to fluctuations, the 
formation of a block is inevitable. As (3 is increased beyond the transition line into 
the Shock-Low phase, the possibility rises that this block will be stabilized by the non- 
conserving dynamics. A spontaneously created block of one species, which now has a 
stable macroscopic size, hinders particles of the other species and breaks the symmetry 
between the two species. 

As (3 is increased from its value at the boundary line between the symmetric shock 
phase and the Shock-Low phase, the size of the block of the majority species increases. 
This is due to the coupling between the ejection rate j3 and the effective injection rates, 
which at this region of phase space serves to increase a + (in mean-field this can be seen 
from Eqs. ITT ^ 125 )1 . At some value of (3, the block reaches the size of the entire system, 
and there it stays for some range of (3. This, in fact, is the High-Low phase of the 
system. As (3 is increased further the size of the majority block shrinks back, and the 
system is again in the Shock-Low phase. The transition from the Shock-Low phase to 
the High-Low phase at some (3, and the re-entrance to the Shock-Low phase at some 
higher (3, occur at these points where the size of the majority blockage becomes identical 
with the size of the system. 

Towards the boundary line between the Shock-Low phase and the low-density phase 
the size of the majority block vanishes. The existence of a reminiscent block which yields 
the asymmetric Low phase, as discussed in jH], can be either attributed to fluctuations 
of the localized shock, or to an alternative mechanism. 
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To make the blockage picture more quantitative, we calculate the size of the 
majority block in each phase. This is done within mean-field for the bulk-decoupled 
case. For simplicity we take up = 0, where the profiles of the blocks are linear. The 
block size is then just 1 — x s , where the shock position x s is given in (|24jl . In figure 
EH (a) we plot the size of the block as a function of (3 at constant Q and a. Using the 
picture described above, one can identify the phase boundaries. The value of (3 at which 
the blockage first disappears corresponds to the Shock-symmetric to Shock-Low line. 
The two values between which the blockage spans the system are identified as the two 
lines between the Shock-Low and High-Low phases. Finally, the higher (3 at which the 
blockage disappears completely corresponds to a transition into a low-density phase. It 
is easy to verify, by comparing with the mean-field phase diagram, that these are indeed 
that transition points between the phases. Note that the Low asymmetric phase escapes 
this picture. 

It is also possible to describe in terms of the block size (or alternatively the shock 
position) the fact that asymmetric phases disappear as Q is increased. The non- 
conserving dynamics in the bulk of the system serves to sustain the localized shock. 
Keeping a constant, for example, the increase in the amplitude of the non-conserving 
rates drives the position of the shock out of the system, thus decreasing the maximal 
size of the blockage. In terms of phases, this would decrease the segment on the f3 axis 
in which the shock is localized at the 'entry' boundary {i.e. the High-Low phase), down 
to a point where the shock cannot get so far and the phase disappears. Beyond this 
point, as depicted in figure El (b), the position of the shock is driven back towards the 
'exit' boundary, thus reducing the Shock-Low phase until it is finally gone. 

8. Summary 

In this work a two-species one-dimensional model, with dynamics which is not conserving 
both at the boundaries and in the bulk, has been studied. The dynamics is symmetric 
under charge exchange and left-right reflection. By definition the non-conserving 
dynamics in the bulk of the system acts to diminish the difference between the densities 
of the two species. Nevertheless we have found that the symmetry between the two 
species can be broken even in the presence of bulk non-conserving processes. 

The mean-field phase diagram, obtained for the case where the bulk dynamics of the 
two species becomes decoupled, exhibits three phases in which the symmetry between 
the species is broken. One of these is unique to the case where the bulk dynamics 
is not conserving, and results from a localization of shocks in the density profiles. All 
asymmetric phases reside in regions of phase space where symmetric low-density profiles 
are another fixed point of the mean-field dynamics. However, stability arguments shows 
that it is the asymmetric solution which should be expected to survive fluctuations. 
Indeed, comparing with Monte-Carlo simulations, two asymmetric phases are confirmed. 
The third, in which the average density of both species is below | is more difficult to 
determine. 



SSB in a non- conserving model 



22 



In contrast to the bulk- conserving case, in this model the density profiles are 
generally not flat. In particular, localized shocks may be generated in the bulk of 
the system. In the general case, when the particle current of one species depends on 
the density of the other, a localized shock in the density profile of one species induces a 
shock in the other. 

In the asymmetric phases, as well as in the shock symmetric phase, typical 
configurations can be described in terms three blocks. The leftmost block has a high 
density of negative particles, the middle block has a low density of particles of both 
species, and the right block is mainly occupied by positive particles. This observation 
serves to define a toy model which describes the dynamics of the system in terms of a 
two-dimensional random walker. Solving the toy-model yields an exact solution for the 
case where the exit rates and the non- conserving rates are taken to zero. The results 
coincide with the ones obtained in mean- field at this limit. For the general case a more 
qualitative picture emerges, which serves to describe the phase transitions in the model 
in terms of the block sizes. 

The bulk of this work, as well as of those which studied the Bridge model, has 
focused on the case in which the dynamics of the two species is decoupled in the bulk. 
The other, more general case, was studied only by numerical means, both on the mean- 
field level and in Monte-Carlo simulations. This enabled us to observe induced shocks. 
A more detailed study of this case by analytical means should shed more light on this 
phenomenon. 

The mean-field phase diagram of the NC-TASEP is expected to be exact [13J, while 
Monte-Carlo results suggest that this is not the case here. Such is also the case in the 
Bridge model. It should be interesting to study the correlations which build up in the 
system, taking it away from the mean-field description. 
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Appendix A. Mean-Field analysis of the case uj d > 

In this Appendix we discuss the construction of the phase diagram for the bulk-decoupled 
case q = 1,ua = wx with ujd ^ 0. This case corresponds to the NC-TASEP with non- 
equal attachment and detachment rates (the case K ^ 1 in fIJEEE])- The P nase 
diagram of the model in the case ujj ^ exhibits two symmetric phases, symmetric 
low-density phase and shock symmetric one, and three asymmetric ones, High-Low 
density phase, Shock-Low phase, and low-density asymmetric phase. The profiles do 
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not have constant slopes as in the up = case. Starting from left boundary density a, 
the density at the right boundary resulting from the left characteristic is given by 



2-uW ( ( 2 - 4a ~ 2cm ) cxn ( 4«+2q M +^(2+») i -2 
W ( \ \ u V \ u J J 

pj(a) = ; ^- - J -. A.l 

HNK ' 2(2 + u) V ' 

The respective expression for the density at the left boundary resulting from the 
right characteristics starting from density 1 — /3 reads 

PV(P) = ^ "\ 11 ■ (A.2) 

2(2 + u) 

Our aim is to use the known phase diagram of the NC-TASEP [11 to construct the 
phase diagram of the two-species model, as it was done in section 13.11 To this end, let 
us define the two transition lines in this phase diagram. The first is the transition line 
between the high density phase and the localized shock phase, 

2(1 + u ) + uW ( (4a - 2 - 2M+2aM) exp ( - 4a+2+2 "- 2 - +n(2+M)2 ) ) 

Ma) = ^ \ , V LL . (A.3) 

2(2 + u) 

where Wk(z) is the Lambert- W function. This equation defines the line only for a < 1/2. 
The second transition line separates the shock phase and the low density phase. This 
line is given, for a < 1/2, by 

2 _ uWi ( (2-^-2^) / _ 4 q+2 a u+n(2 + u)l -2 \ \ 

AlW = ' l - "\ 11 (A.4) 

2(2 + u) 

Both boundary lines are continued for a > 1/2 by Phs(^) an d /3sl(|), respectively. 

Let us recall the procedure in which one obtains the phase diagram for the two- 
species model. For each phase, one obtains from (J13)) the effective boundary rates, a + 
and a". This requires knowledge of the four boundary currents, jt = Pb(^ ~ Pb) an d 
j b = m b {l — m b ), where b = 0, N for the left and right boundaries. The boundary lines 
are then obtained from comparing the effective boundary rates with the corresponding 
transition lines of the NC-TASEP. 

We do not repeat the analysis here in such details as it was done for the case uoo = 0. 
Instead, we give for each phase the four boundary densities, needed to calculate the 
boundary currents and the effective rates. In addition the conditions on the effective 
rates, which define the phase boundaries, are given in terms of /?hs and /?sl of Eqs. 
(IA.3j) and (|A.4J) . The asymmetric low-density phase is omitted, as we could not obtain 
the boundary densities in this phase in a closed form. The boundary lines for this case 
were obtained numerically. 

Low density symmetric phase. Here, po = tun = ct + , Pn = mi = P^n{o- + ) as in 
flA.lj) . The condition for the existence of this phase is 

(3 < (3 Sh (a + ) . (A.5) 
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Shock symmetric phase. The boundary densities are given by p = m N = a + , 
Pn = m i = 1-/9. The conditions for this phases existence are 

MO </3<M« + ) • (A.6) 

Shock - Low asymmetric phase. Let the positive particles be in the shock phase. 
Then p = a + , p^ = 1 — (3. The negative particles are in the Low phase, where = a" 
and m\ = p^(a~) as given by equation lA.il The phase exists in a region in phase space 
where 

MO < (3 < MO • (A.7) 

High - Low asymmetric phase. As before it is assumed that the positive particles 
are in the majority phase: p = Pi\(3) according to IA. 21 and p^ = 1 — (3. The negative 
particles are in the Low phase: = or and m\ = p^(a^). The High-Low phase 
exists where 

MO < (3 < MO • ( A -8) 
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